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NOMENCLATURE 


A k 


homogeneous polynomial of Kth degree in Z, S 


A s 


speed of sound where air speed equals zero 


A 0> A l> 
A 2 > A 3 


parameters of polynomial PI 


a 

C 


D, E 
h 

I(2.,m,n) 
i,j ,H , m 
k 


M 

N 


n 

P 

PI 

PN 

q 

s 


s,z 


T (n) 

V 

X,Y 

Z 


speed of sound at arbitrary point in the flow 

parameter of polynomial PI 

infinite series of multiple integrals 

function of the gas constant k 

symbolic coefficients of integral of quantity E 

summation indexes 

ratio of specific heats 

Mach number 

degree of polynomial PN 
iteration number 
polynomial in Z, S 
polynomial to be integrated 
polynomial of degree N 
speed of sound 
numerical value of S 
lower limits of integration 
upper limits of integration 
nth integral of quantity E 
velocity vector 

Cartesian coordinates in physical plane 
numerical value for Z 


CD 

Cartesian 

CD 

O 

V* 

O 

constants 

e v 

numerical 


numerical 

VI, v,? 

indexes 


potential 


coordinates in pseudologarithmic plane 
in pseudologarithmic plane 
value for 6 
value for X 

function of a compressible fluid flow 


IV 



COMPUTATION OF THE POTENTIAL FUNCTION OF A 


COMPRESSIBLE FLUID NEAR MACH NUMBER 1 
Howard Tashjian 
Ames Research Center 


SUMMARY 


A procedure is described for finding the potential function of a 
compressible fluid flow according to Bergman f s method of orthogonal functions. 
In Bergman's formula for the potential function an infinite series occurs, the 
terms of which are multiple integrals. These multiple integrals have a singu- 
larity at Mach number 1. The problem arises of computing the multiple 
integrals as the Mach number approaches 1. 

To calculate the integrals, it is necessary to use a special 
computational system. A method is described whereby a digital computer 
equipped with a Formula Manipulation Compiler known as Formac can be used to 
compute and evaluate the multiple integrals. 


INTRODUCTION 


In reference 1 a method is described for computing flows of 
compressible fluids. Research undertaken with the object of computing flow 
patterns of compressible fluids has led to large expressions to be manipulated 
symbolically. The expressions that have arisen in the formula for the poten- 
tial function are so large that a human could not analyze them mathematically 
(ref. 1, p. 305). In the present investigation, attention has been directed 
to the use of a digital computer for performing the mathematical analysis 
(ref. 2) . 

The extensive storage capacities of the digital computer, combined with 
the development of a new programming language known as Formac, enable it to be 
used to perform the mathematical analysis of problems of the type under con- 
sideration (ref. 3) . As a consequence of the algorithm for performing symbolic 
integration of a power series, the operations involved in accomplishing the 
symbolic integrations of algebraic terms arising in the formula for the poten- 
tial function of a compressible fluid can be translated into a Formac program. 
The aspect of this problem which makes it so attractive for computer solution 
is that the mathematical analysis can be translated into a computer program. 

The purpose of this report is to show how the computer can be used to 
perform the mathematical analysis involved in the evaluation of the formula 
for the potential function of a compressible fluid as the Mach number 
approaches 1. With this end result in mind, a Formac program was written to 
perform symbolic double integrations of polynomials of arbitrary degree 
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combined with the facility for performing an arbitrary number of iterations. 
To compute the potential function for different Mach numbers, the user need 
only supply the computer with the appropriate polynomial, and specify the 
desired number of iterations. 


ANALYSIS 


For two dimensions, the potential function <|> of compressible fluids 
when considered in the physical plane satisfies the nonlinear partial 
differential equation 


3 2 4 

i - 

)<P\ 2 

2 34 34 3 2 4 

+ ifi 

i -±(; 
a 2 V 

\ 2 1 

)4\ 

3X 2 

L a 2 

J 

2 3X 3Y 3X 3Y 
a^ 

3Y 2 

*Y/ J 


= 0 


CD 


where 


2 = a 2 


a^ = a 


k - 1 


'/ ,2 , >2 
/ jhM /8<M 

\3X/ + \3Y,/ 


(ref. 4, pp. 239, 241) 


( 2 ) 


and a s and k are constants. 


Equation (1) is first linearized by transformation to a pseudologarithmic 
plane with Cartesian coordinates 0 and 


A = 2" log 


1 - (1 - Mp 


211/2 


1 + (1 - M 2 ) 1 / 2 


1 + h Cl - M 2 ) 1 / 2 


1 - h(l - M 2 ) 1 / 2 


1/h 


(3) 


The coordinate 0 is the angle formed by the velocity vector 


u i0 
V = qe 


(4) 


and the positive X axis, and q is the speed (ref. 1, p. 302). 

. 1/2 


-(H 4 ) 


C5) 


and 


M = 


q 

(k - 1) q 2 / 2] x / 2 


( 6 ; 


where M is the Mach number (ref. 1, p. 302). The physical and 
pseudologarithmic planes are compared in figures 1 and 2. 
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The following formulas are also 
obtained : 



Figure 1.- A flow in a channel. 

X 




( 7 ) 



q 2 ctn 2 e 
1 + ctn 2 0 

q 2 

1 + ctn 2 6 


( 8 ) 


(ref. 1, p. 304) 

(9) 


For flows satisfying equations 
(1) through (9) , the potential function 
has a singularity on the sonic line. 

In reference 1 (p. 305), an analytic 
expression for this type of singularity 
has been determined in the pseudologa- 
rithmic plane. In this plane, the 
formula for the potential function is 

4> = | [log (Z - Z 0 ) + log (S - S 0 )] + D 

( 10 ) 


Figure 2.- The image of the flow in figure 1 in 
the pseudologarithmi c plane. 


where 


Z = 0 + iX 


S = 0 - iX 

Z 0 = 0q + 


ref. 1, p. 305 


S 0 = 0q - / 


( 11 ) 


For the quantity E is obtained an infinite series whose terms are multiple 
double integrals. In the following section is presented the computational 
procedure for the formal integration and numerical evaluation of these 
multiple integrals. 
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Computational Procedure 

For the quantity E appearing in equation (10) are obtained multiple 
integrals of the form 


T 


( 1 ) 


T 


( 2 ) 


T 


(3) 



PN(Z,S)dZ dS 


T ( - 1 - ) PN(Z,S)dZ dS 


T(2) PN (Z , S) dZ dS 


-\ 



T(n) 




T (n-l)P N (Z,S)dZ dS 


(ref. 1, p. 305) 


The subscript of T in equations (12) denotes the iteration number 
referred to later. An expression similar to equations (12) is obtained for 
the quantity D in equation (10) (ref. 1, p. 305). In this report is dis- 
cussed the formal symbolic integration and evaluation of equations (12) . The 
results presented were obtained by using the IBM 360/50 computer and the 
Formac programming language. The Formac programming language is a practical 
system for doing nonnumeric mathematics on the computer. 


In reference 5 a method is given for determining the polynomial PN. If 
PN is a polynomial of degree N in Z,S, then, as formal computation shows, 
T( n ) is also a polynomial in Z,S; namely. 


T 


(n) = 


E 


I (&,m,n)Z^S m 


(13) 


where I(£,m,n) are symbolic coefficients dependent on Zq and Sq and i and 
m are the powers of Z and S, respectively, and n is the iteration number. 
In equation (13) the products Z^S m are ordered in such a way that at first 
those terms appear for which £ + m = p = maximum. Then the process is 
repeated for terms of degree (p - 1), (p - 2) , etc. According to this 
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procedure, T can be expressed as a sum 


T ( n ) = A p + A p-i + • • • + A o 


(14) 


where is a homogeneous polynomial of the Kth degree in Z,S. Then in 
every the terms are ordered in descending powers of S. 

A recursion procedure for determining the coefficients I(£,m,n+1) from 
the known coefficient I(£,m,n) will be discussed next. From the definition 
of equations (12), it follows that 


T (n+i) 


■ r r « « 

Z 0 •'so 


(15) 


The first step in determining I(£,m,n+1) is to form the polynomial 


T 


(n)PN 


^ ^ (n) ^ ,m,n)Z £ S m 

£, m 


EiCiJ.DzV 

i »j 


where PN has been denoted by 


(16) 


PN 


X\(i,u)zV 

i» j 


(17) 


The polynomial represented by equation (16) can be determined by the use of 
Formac, and the result written in the form 


T(n)PN 



(£,m,n)Z £ S m 


£,m 


(18) 


The next step in determining I(£,m,n+1) is to perform a double integration on 
the right member of equation (18) . 


' (n+i ) 


■/7‘E 

Z 0 Sq £,m 


P(£,m,n)Z £ S m 


dZ dS 


( 19 ) 


Formac can be used to perform this double integration if the right side of 
equation (19) is replaced by 



(n+l) 


£ 

£,ni 


gCj-AB) 


(£ + 1) (m + 1) 


z* + 1 s m+1 


„£+l c m+i 
Z 0 b 


£+1 m+1 _£+i m+1 

- L oq + z 0 t>0 


( 20 ) 


and the terms are arranged in descending powers of 


,£+l 


s m+1 


Computer Results 


If a set of values Zy = 0 V + iA^, Sy = 0 V - iA^, y = 1, 2, . . . , p is 
given, the values of T( n ) at these points can be determined by the use of 
Formac. The integrated expressions, as well as the expressions to be inte- 
grated, become very long. Therefore, only a few typical integrals and inte- 
grands will be presented, together with some numerical values of the integrated 
expressions . 


In equations (21) a third degree polynomial is shown to be integrated. 
The name of the polynomial is PI . 


P 1 ( 3 » 0 » 

1) 

= 1/8 

A3 

Pl(2,li 

1) 

= 3/8 

A3 


Pl(2,0,l) 


1/4 A2 + 

3/8 A3 

C 

PL ( 1, 2, 1> 

= 

3/8 A3 



Pl(l, 1,1) 

- 

t/2 A2 + 

3/4 A3 

C 

“P 1 x~n 0 , FT - 


17Z“AT + 

“1/2 A 2 

2 

378 A3 C 

P 1 ( 0 , 3 , 1 ) 

= 

1/8 A3 



PKO, 2, 1) 

= - 

1/4 A2 + 

3/8 A3 

C 

Pl(0rl, 1) 

= 

1/2 A 1 * 

1/2 A2 

2 

C + 3/8 A3 C 

P1(0, 0, 1) 

= 

AO + 1/2 

C A1 + 

2 

1/4 A2 C + 1/8 


} ( 21 ) 


The first subscript of PI indicates the power of Z, the second 
subscript the power of S, and the third subscript the iteration number. 

In equations (22) are shown the integrals of equations (21) . The name of 
the integral is I. The first subscript of I indicates the power of Z, the 
second subscript the power of S, and the third subscript the iteration 
number. 
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( 22 ) 


For the numerical evaluation of equations (22) the parameters of the 
polynomial PI are assigned the following values: 

~\ 

A 0 = 37290.0 

k\ = -0. 8426261. 10 7 


A 2 = -0. 56560455. 10 10 
A 3 = -0.33684736. 10 13 
C = 0.0045 


(23) 


Z 0 = -0.00225 
S 0 = -0.00225 


J 
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For convenience, the formulas for Z and S are reproduced. 

Z y = e v + iA c (24) 

S p = 0 V - iAg (25) 

Five distinct values are assigned to 0 V as represented by formula (26) . 

e v = 0.1(v - 1) , 1 < v <- 5 (26) 

For each of the five values of 6 V , four values are assigned to A^ as 
represented by formula (27) . 


A = -0.003 + 0.0005(C - 1) , 1 <. ? <. 4 


(27) 


Twenty values of Zy,Sy are obtained, yielding twenty values for equations 

( 22 ). 


The numerical values of equations (22) are shown in equations (28) . The 
names of the numerical values are denoted as SUM(v,5,l), where v,5 are as 
defined above, and the third subscript 1 denotes the first iteration. Each 
of the SUM (v, 5,1) in equations (28) represents a numerical value of equations 
(22) for a given 0v,Ag. To the right of the first equal signs in equations 
(28) are shown the numerical values of equations (22) . A second equal sign 
appears in equations (28) beginning with the f j fth SUM. This second equal 
sign, together with the symbol I which follows it, represent /^l. There- 
fore, the first parts of the numerical values are complex, and the second 
parts real. 
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SUM ( l r 1 r 1 ) 


02238094 


soMtt. 2 , r r- 

— 

.0023^062 



SUM<1,3,1> 

= 

.00225773 



SUM( 1,4,1) 

- 

.01852636 



SUM<2,1,1) 

= 

- .13633494E-C9 =1 

+ 

8175.47963 

SUM(2,2,1> 

= 

- .26372815E-09 =1 

+ 

18686.2335 

SUM ( 2,3,1) 

= 

- • 64614535E-10 =1 

* 

29191.5161 

SUM! 2(4,1 I 


- . 11833733E-09 =1 


39688.1661 

SUMOflf 11 

= 

- • 1542331 5E-C8 =1 

+ 

126027.486 

SUM ( 3 , 2 , 1 ) 


- .41026835E-08 =1 

+ 

294388.321 

SUM 13,3,1) 

= 

. 69750 1 61E-09 =1 +■ 

462727.212 

SUMO, 4,1) 

= 

. 10666321E-08 = 1 + 

631031.523 

SUM { 4 r 1 » 1 > 

= 

- • 72820669E-08 =1 

+ 

633532.947 

SUWt 4 v 2 y 1 ) 

= 

- • 1 035 8456E-C8 =1 


1486036.45 

SUM ( 4 , 3,1) 

= 

. 31876594E-08 =1 + 

2338490.55 

SUM( 4,4,1) 

= 

- .68060566E-C8 = 1 


3190866.83 

5UK( 5 fit IT 

= 

'.'42211 8IE- 07 =T + 1997320.03 

SUM ( 5 ? 2 » 1) 


- .51 7938 13E-07 =1 

+ 

4691847.6 

SUM<5,3, 1 ) 

= 

- . 78278551F-C7 =1 

+ 

7386287.32 

S UM { 5 r 4 f 1 ) 

= 

- .54500574E-07 =1 

+ 

1.0C805886E*07 


(28) 


The polynomial for the second iteration is obtained by multiplying 
equations (22) by equations (21). The multiplied expressions become so large 
that only two typical terms will be shown in equations (29) . 
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(29) 

The integral of equations (29) for the second iteration becomes so large 
that only one term will be shown by equation (30). 
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CONCLUSION 


The availability of the large-scale digital computer, together with the 
development of symbolic programming languages such as Formac, permits the use 
of the mathematical solution to problems in compressible fluid flow theory. 

It has been shown that the memory capacity of the computer can accommodate 
many thousands of symbolic terms arising in the formula for the potential func- 
tion of a compressible fluid flow. Furthermore, the possibility of numerically 
evaluating these many thousands of terms has been demonstrated. 

The results presented in this report suggest that the computational 
techniques employed here permit the solution to problems in compressible fluid 
flow theory which could not be accomplished by other means. 
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